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1 Introduction 



Rare decays of the B meson proceed via a variety of diagrams. They provide a unique 
opportunity to study the Kobayashi-Maskawa matrix elements, \Vcb\, \Vub\, \Vts\ and \ Vtd\, 
as well as their relative phases. In addition, decays involving loop diagrams are interesting 
because of their sensitivity to effects of new physics. Such rare B meson decays can be 
studied in detail at the recently commissioned e^e~ -B-factories, KEKB and PEP-II. 

By the end of December 2000, the Belle experiment at the KEKB accelerator |^ 
had accumulated 10.4 fb~^ of data on the T(4S) resonance, corresponding to 11.1 x 10^ 
BB pairs. Based on this data set, results for the following rare decays are briefly re- 
viewed in this paper; 1) Cabibbo suppressed decays B D^*^K^*\ 2) charmless two- 
body decays B tttt, Kn, KK, 3) charmless decays involving the rj' meson B 
ri'h(h = K^,Tc^, K^)^ and 4) a search for flavor changing neutral current (FCNC) de- 
cays B K^*H+e-. Resuhs for additional rare decay modes can be found in other 
contributions to this conference 0]. 

In all the analyses, B candidates are reconstructed using the beam constrained mass, 
^bc = y^-E^eam " p|> ^ud the energy difference, AE = Eb - -Ebeam- Here, -Ebeam = 
^/s/2 ~ 5.290 GeV, and pb and Eb are the momentum and energy of the reconstructed 
B in the T(4S) rest frame, respectively 0. Normally we compute AE assuming a vr 
mass for each charged particle from B decays. This shifts AE by about —50 MeV for 
each meson. Other common analysis techniques, the particle identification (PID) 
cut and the qq background suppression cut, are described in appendices. Throughout 

this paper, the inclusion of charge conjugate states is implied, except for the direct CP 
°In the analysis for the Cabibbo suppressed decays, described in Sec. the beam constrained mass is 
calculated in the laboratory frame as nibc — \J {E^g^Y ~ IPb'''P- Here p^g^ is the B candidate's laboratory 
momentum vector and i?^'' — -^{s/2 + Pee ■ Pb'^), where s is square of the center of mass energy, and 
Pee and Eee are the laboratory momentum and energy of the e^e~ system, respectively. 
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violation measurement described in Sec. |^. The first and second errors in results represent 
statistical and systematic errors, respectively. 



2 Cabibbo Suppressed Decays 

Direct CP- violating asymmetries due to interference between b ^ c and h ^ u transition 
amplitudes in the Cabibbo suppressed B~ D^K~ decay is a theoretically clean way 
to determine the 03 angle of the unitarity triangle [Q. As a first step of this program, 
we have searched for the Cabibbo suppressed processes D*^K^, D*^K~, 

B^ D^K^ as well as B^ —>■ D^K~ . In the tree level approximation, their branching 
fractions are related to those of the Cabibbo favored B D^*\^ counterparts by 

R n(*) T{-\ 

R - ^ ^W,-) - tan^ ecifK/Uf - 0.074, (1) 

using the Cabibbo angle Oc and the meson decay constants fK{-K) 0- The only Cabibbo 
suppressed decay observed to date is B^ D^K^, reported by CLEO p with R = 
0.055 ±0.014 ±0.005. 

In the analysis, D mesons are reconstructed using the decay modes, D° — > K~7r~^, 
K-TT+n^, K-n+n+n- and D+ K-TT+n+, K°s^+, Kln+n+n-, K~K+n+. For and 
D*^ reconstruction, D*^ D^tt^ and D*^ —>■ D^n^, D^n^ decays are used. If multiple 
entries are found in one event, the best candidate is selected based on a determined from 
the differences between measured and nominal masses of rriD, rribc and, when appropriate 
m£)* — rriD and m^ro (vr^ from the D* decays). To suppress the much more abundant 
B D^*^!!^ decays, a relatively tight PID cut IZk > 0.8, which gives 76.5% efficiency 
for K tracks and 2% n —>■ K fake rate, is applied for the prompt hadron track from the 
B decay (see Appendix A). 

^ The relation assumes the vahdity of factorization and flavor-5'C/(3) symmetry. 
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Figure shows the AE distributions for all B D^*^n- {TZr < 0.8) and B D^*^K- 
{TZk > 0.8) enriched samples. In all the B D^*^K channels, signal peaks are clearly 
seen at AE = — 49MeV. Table |l] summarizes the results. The R ratios, defined in Eq.(|l]), 
are calculated from the extracted yields of 5 ^ D^*'^'rr^ and B D^*^K^ decays by 
taking into account the difference in the reconstruction efficiencies. In all cases, the R 
ratios are consistent within errors with the expectation in Eq.(|lD. These are the first 
observations of the B D^K^ , D*^K~ and D*^K^ decay processes |^ 

We have also searched for the B^ D^K*^ decay mode. The K* mesons are recon- 
structed using the clean Kg7T~ channel. The first observed evidence of this decay mode is 
shown in Fig. |^. The signal yield is 15.0 ± 4.6 with a statistical significance of 4.6 a. The 
branching fraction normalized to that of B^ D^n^ is found to be 0.1 16 ±0.036 ±0.015 
(preliminary) Q. 

3 B ^ TTTT, Ktt, KK Decays 

Branching fraction measurements of the B ^ mT, Kn and KK decays are an important 
first step toward indirect and direct CP violation studies with the tth and Kn modes, 
which are related to the angles 02 and 03 of the unitarity triangle, respectively. Moreover, 
assuming isospin and SU{3) symmetry, recent theoretical work suggests that branching 
fraction of these modes can be used to extract or place limits on the angles 02 and 03 ^ . 

We have analyzed the data for the B^ tt+tt", K^tv^, K^K~, K^tt^ and 5+ tt+tt^, 
K^n^, K^n^ decays. The Belle detector is equipped with a high momentum PID system, 
which gives clear separation between the charged vr and K mesons in these final states. 
The 71 and K mesons are distinguished by the PID cut 7l-n-(K) > 0.6, which gives 92.4% 
Results on B ^ D'^*^K modes were updated after the conference. 

The ratio similar to Eq.(0), R = B{B- D"K*-)/B{B- D°p-), is not presented here because 
our analysis for the Cabibbo favored counterpart B^ has not been completed yet. 
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efficiency and 4.3% fake rate (true vr fakes K) for vr mesons and 84.9% efficiency and 10.4% 
fake rate (true K fakes vr) for K mesons, respectively (see Appendix A). The dominant 
background from the continuum qq process is suppressed using the hkehhood ratio method 
described in Appendix B. 

Figure ^ shows rrifyc and AE distributions for the 7t~^tt~, K~^7i~ and Kg7T~^ modes. 
Each distribution is ffited to a Gaussian signal on top of a background function, which 
is modeled by the ARGUS background function for rribc and by a linear function for 
AE. The AE ffi results are used to determine the signal yields. For the vr+vr^ and K^n^ 
modes, the cross-talk from each mode is not negligible, and thus the ffis treat the size of 
both components as free parameters. The cross talk found in these fits are consistent with 
expectations. Figure § shows the mbc and AE projections for the tt^tt^, K^tt^ and K^tt^ 
modes. For these modes, since the AE distribution has a long tail, a two-dimensional fit 
is applied to the mbc and AE distributions. For the 7r"^7r° mode, since the cross talk from 
K^n'^ is significant and separated by less than la in AE, the ffi includes the expected 
K'^tt'^ component. 

Results of the branching fraction measurements are summarized in Table |^. The 
charge averaged branching fractions for B vr+vr", K^tt", i^+vr^, K^tt^, and K'^tt^ 
are measured with statistically significant signals. Our results are consistent with other 
measurements con&m that B{B^ —>■ Ji+vr^) is larger than B{B^ tt+tt^), and 

indicate that B{B^ h^n^) and B{B^ K^tt^) is larger than expected in relation to 
the B^ h^TT~ and B^ K^tt^ modes based on isospin or penguin dominance 0. 

We have also searched for direct CP violation in B ^ K'^tt^ and B K^ir^ modes. 
With the high momentum PID at Belle, the dilution of the CP asymmetry due to double 
misidentification is as small as ~ 0.5%. Inherent asymmetries in the Belle detector are 
determined to be less than 2% based on the yield difference between — > K^tt^' and 
Results shown in Table were updated after the conference. 
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D° — s> K^TT^ decays. Figure ^a) and b) shows the AE distributions for the K~it~^ and 
K^Ti^ final states. The signal yields, extracted with the same procedure as the branching 
fraction measurement, are 27.7^6 i for K^n^ and 25.4^63 for K~^7i~, giving 

a partial rate asymmetry Q AcpiK'^Tv^) = 0.043 ± 0.175 ± 0.021 and its 90% confidence 
level interval —0.26 < AcpiK^^n^) < 0.35. Figure ^c) and d) show the same distributions 
for the K-7r° and K+n^ final states. The yields are 18.3^^^ for B' K~7t° and 17.61^1 
for B^ — > K'^TT^, giving AcpiK'^jr^) = O.OlOlg^g^ (only statistical error is shown for this 
mode). It should be noted that these results on Acp are preliminary. 



4 B ^ rjh Decays 



Measurements by CLEO [10| and BABAR indicate that the branching fraction for B — >■ 
7] K is significantly larger than theoretical expectations. Here, we present preliminary 
results for the B^ V ^ 'q'rc'^ and B^ r]' decay modes. The r]' mesons are 
reconstructed using rj' rj-R-R, ^ 77 and rj' —> p7 decays. The qq background is 
suppressed by a likelihood ratio cut similar to that used for the B tttt. Kit, KK 
analysis. The likelihood contains SFW and cos 9b (see Appendix B). The prompt charged 
K/n from the B decay is separated by the PID cut TZ^^i^K) > 0.6. 

Figure ^ shows the rrifyc and AE distributions for the signal candidates. By com- 
bining results for the two sub-decay modes of the rj' meson, the observed signals have 
statistical significances of 10.9cr and 5.9cr for the t]' K'^ and r]' K'^ modes, respectively. 
Their branching fractions are found to be B{B^ — * rj' K^) = (6.8 li;! -o'l) x 10~^ and 
B{B° ri'K'^) = (6.4 x 10~^ In the rj'ir-^ mode, no statistically significant 

signal is seen, and we set an upper limit of B{B^ r] -n^) < 1.2 x 10^^ (90% C.L.). 
These results are consistent with the previous measurements. 



^ Here the partial rate asymmetry Acp is defined as, Acp = j^|gZ^jj_|_j^|gZ^j| ■ 
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5 Search for B K^*k^£- Decays 

The FCNC decay, B K^*H+i-, proceeds via loop dia Eframs, and therefore is sensitive 
to new physics such as charged Higgs and SUSY. We have searched for the FCNC decays, 
^0 ^ K*H+i-, B+ 5° ^ K^i^i-, B+ K+i+i-, in both electron and 

muon channels. 

The mesons are reconstructed using Kg tt+tt^ decays. The K* mesons are 
reconstructed using K*^ K'^it~, Kgir^ and K*^ — > K^tt^, K^tt^ decays. Electrons 
with piab > 0.5 GeV/c and muons with piab > 1-0 GeV/c are selected with pion fake rates 
expected at a level of 0.3% and 1.7%, respectively. 

Background from qq processes is suppressed by a likelihood ratio cut (see Appendix B) , 
where the likelihood contains cos 6b, cos 6 da and a Fisher discriminant constructed from 



virtual calorimeter variables ||TT[ and the normalized second order Fox- Wolfram momen- 
tum i?2- The largest source of BB background are events where both B mesons decay 
into iuX. This source is suppressed by cutting on a likelihood ratio based on visible 
energy and cos 6b- The large background from J/ijj ^ i~^£~ and ip' — > decays is 

vetoed by applying the di-lepton mass cuts, —0.15 < Mee — < 0.07GeV/c^ and 

—0.10 < M^^ — Mjj^^^i-^ < 0.05GeV/c^. The di-electron mass is required to be larger 
than 0.1 GeV/c^ to avoid background from -k^ Dalitz decays and 7 conversions. 

No statistically significant signals have been observed yet, and upper limits at 90% 
confidence level are set for each mode, which are shown in Fig. |^ along with a comparison 



to theoretical expectations [0 and the results from previous experiments ||r3|'||14||'|[T5 



Our upper limits are more restrictive than those of previous experiments, except for the 
K*^fifi channel. 



7 



6 Summary 



In this paper, we briefly review results for some B meson rare decays from the Belle exper- 
iment, based on 10.4 fb~^ data collected on the T(4S) resonance. They are summarized 
as follows, 

• Cabibbo suppressed decays B D^*^K have been observed with nearly the expected 
size relative to B ^ D^*^7r decays. We claim the first observation of the four modes, 
B D+K-, D*^K-, D*+K- and D^K*-. 

• Our results for the B —>■ tttt, Kit, KK modes are consistent with other measure- 
ments, confirm the small tttt/Kt: ratio, and indicate that the K^t:'^ mode is larger 
than theoretically expected. 

• The large branching fractions for the B rj' and rj' modes are also confirmed. 

• Improved upper limits on the FCNC decays B — > K^*^i~^i~ are obtained except for 
the K*^ijl~^ijl~ channel. 

In the near future, we anticipate results on various rare B meson decays with much higher 
statistics and better sensitivity than ever before achieved. 
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Appendix A: particle identification cut 

A clear separation of charged tt and K mesons is essential to find some rare decay signals. 
We combine outputs from three sub-detectors, specific ionization loss in the central drift 
chamber (dE/dx), time-of-fiight measured by scintillation counter arrays (TOF) and light 
yields measured by the aerogel Cherenkov counter arrays (ACC), into a likelihood for each 
particle assumption CK{n) = I^kI-k) ^ ^^^) ^ ^"kI^Y "^^^ particle identification is then 
performed by cutting on the likelihood ratio TZk{-k) = I^k{-k)I{C.k + C.-k)- For prompt tt 
and K mesons from B decays, TOF does not provide useful separation, and only dE/dx 
and ACC information are used. The PID cut efficiency and fake rate are measured using 
kinematically selected D*^ — > D^tt"*", — > K'ti^ decays. 



Appendix B: continuum suppression cut 

To distinguish BB events from the qq background. Belle has developed a new event shape 
variable called "Super Fox- Wolfram" . The usual Fox- Wolfram moments are defined as 
Hi = Sij |Pi||p^|-Pi(cos6'ij), where the indices i and j run over all final state particles. 
Pi and pj are the momentum vectors of the particle i and j. Pi is the l-ih. Legendre 
polynomial, and Oij is the angle between the two particles. In the new method, the 
normalized Fox- Wolfram moments, Ri — Hi/Hq, are decomposed into three terms: Ri = 
Rf + + Rf° = {Hi'' + + H;^°)/Ho, where the indices ss, so, and oo indicate 
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respectively that both, one, or neither of the particles comes from a B candidate. These 
are combined into a six term Fisher discriminant (Super Fox- Wolfram) , defined as SFW = 
Ef=i {aiRr + PiRD, where ai and A are Fisher coefficients and 1=2, A for ai and -Rf°. 
The terms Rf and -Rf=i^3 are excluded because they are strongly correlated with rribc and 
AE. 

We then combine different qq suppression variables into a single likelihood, Cs{qq) = 
Y[i^l(qq): where the -C^^^^-) denotes the signal (qq) likelihood of the suppression variable i, 
and select candidate events by cutting on the likelihood ratio TZg = Csl{Cs + Cqq). The 
variables used depend on the decay mode. For example. Fig. ^ illustrates the separation 
of the signal and the qq background in the case oi B ^ tt+tt^, K^tt^ , K^K^ analysis, 
described in Sec. ^, where the likelihood contains the above SFW , the B candidate flight 
direction (cos6'b), and the decay axis direction {cosOda)- By requiring IZs > 0.8, 97% of 
the qq background is removed while 48% of the signal are kept. In the case of -B ^ tt+tt^, 
K'^tt'^, K^tt^ modes, the event sphericity (S) and the angle between the thrust axis of the 
B candidate and that of all the remaining particles {6t) are also used. In these modes, 
since the continuum background is more severe, a loose cut on cos6t is applied first. 
Then SFW is extended to include cos6't and S, and this extended SFW and cos 9b are 
combined into the likelihood. 
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Table 1: Summary of -B — > D^*'>K results. The extracted yields of D^*^K signals 
Ndk-i their statistical significance E, and the obtained ratios of the branching fractions 
R = B(B D^*^K-)/B{B D^*h-), are shown. 

channel Ndk E R 

B- ^ D^h- 138.4 ± 15.5 11.7 0.079 ± 0.009 ± 0.006 

B^^D+h- 33.7 ±7.3 6.1 0.068 ± 0.015 ± 0.007 

B-^D*^h- 32.8 ±7.8 5.8 0.078 ± 0.019 ± 0.009 

B^^D*+h- 36.0 ±7.1 7.6 0.074 ± 0.015 ± 0.006 
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Table 2: Summary of the B tttt, Kit, KK results. The obtained signal yield (iV^), 
statistical significance (E), reconstruction efficiency (e), charge averaged branching frac- 
tion (B) and its 90% confidence level upper limit (U.L.) are shown. In the calculation of 
the production rates of B^B~ and B^^ pairs are assumed to be equal. In the modes 
with mesons, A^^ and e arc quoted for K^g, while B and U.L. are for K^. Submode 
branching fractions for Kg — > tt+tt" and 7r° — > 77 are included in e. 

Mode Ns E e [%] B [xlQ-^] U.L. [xlQ-^] 

SO^TT+TT- n.ltll 3.1 28.1 0.56 ± 0.04 

S+^tt+ttO 10.4 ili 2.7 12.0 0.78 Ig;! l^' 1-34 

B^^K+n- m.?>tf,i 7.8 28.0 1.93 i[J:i 1°;^^ 

34.9+7.6 7 2 19.2 1.63 1°:!^ 

B+^K^n+ 10.3 It^ 3.5 13.5 1.37 tgj^ 1°;^^ 

B^^K^n'' 8.4 11? 3.9 9.4 1.60 1°;^^ 1°;^ 

5° ^ X+X- 0.2 tli - 24.0 - 0.27 

B+ K+W 0.0 - 12.1 - 0.50 
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AE(GeV) AE(GeV) AE (GeV) AE (GeV) 



Figure 1: The AE distributions for the (a) B' D^h', (b) 5° ^ D+h', (c) B' 
D*^h~ and (d) B^ — > D*'^h~ decay channels, in the iribc signal region (5.27 < rubc < 
5.29GeV/c^). The top figures show B D^*^!: control samples with the PID cut TZk < 
0.8, and the bottom figures show B D^*'>K enriched samples with the PID cut TZk > 
0.8. The points with error bars present the data, the curves show the results of fits. 
The open histograms are the sums of background functions scaled to fit the data and the 
hatched histogram indicates the continuum component of the background. 
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Figure 2: Plots which show evidence for the B~ ^ D^K*~ decay; Distributions of (a) 
for the mf,c signal region, (b) rrihc-AE, (c) Kg7r~ invariant mass for events in the 
rribc-AE signal box and (d) rubc for the AE" signal region. The solid lines in (b) indicate 
the rribc and AE signal regions. 
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Figure 3: The rubc (left) and AE (right) distributions, in the signal region of the other 
variable, for B ^ a) 7r"'"7r~, b) K'^n" and c) KgTi^. The fit function and its signal 
component are shown by the solid and dashed curves, respectively. In the 7r"'"7r~ and 
K'^TT' fits, the cross talk component is shown by the dotted curve. 
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Figure 4: The m^c (left) and (right) distributions for a) tt+tt^, b) K^t:'^ and 

c) KgT^^. For the X'^tt", X mass is assumed for the charged particle. The projection of 
the two dimensional fit onto each variable and its signal component are shown by the solid 
and dashed curve, respectively. In the 7r"'"7r° fit, the cross talk from X+tt^ is indicated by 
the dotted curve. 
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Figure 5: AE distribution for each a) 5° K-n+, b) K+n-, c) K'-n^ 

and d) B+ K+n° mode. Fits are similar to those shown in Fig. |^ and Fig. ^ 
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Figure 6: The rUbc and AE distributions for (top) f]' , (middle) v't^^ 

and (bottom) B^ v' inode. The hatched histograms show the contribution from 
rj' — > p7 sub-decay mode. The sohd hues show the fits with a Gaussian signal on top of 
a background function, which is the ARGUS background function for m^c and a linear 
function for AE. The fits for B^ rj'-K^ include the expected B^ rj' component. 
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Mode : UpperLimit 
KVe" :5.91 x10"^ 



KVe": 2.43x10" ^ 



K*Ve" :11.1 x10"^ 



K*Ve": 20.8x10"^ k 



kW: 8.45x10"^ 



KW": 3.96x10"^ ^ 



K'VV: 10.5x10"^ 



KVV: 8.65x10"^ k 
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Figure 7: Upper limits on the branching fractions for each B — * K^*''ii mode (prehmi- 
nary). Our results are compared to theoretical predictions and upper limits from previous 
experiments. 
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Figure 8: Distributions of the qq suppression variables in the case oi B ^ tt^ti , K^ti , 
K+K- analysis; (A) SFW, (B) | cos^sl, (C) | cos^^al and (D) TZ^. Sohd and dotted lines 
show the distributions for the B decay signal and qq background, respectively. 
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